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ABSTRACT 

We study the behaviour of spherical Voids in Lagrangian perturbation theories L(n), 
of which the Zel'dovich approximation is the lowest order solution L(l). We find that 
at early times higher order L(n) give an increasingly accurate picture of Void ex- 
pansion. However at late times particle trajectories in L(2) begin to turnaround and 
converge leading to the contraction of a Void, a sign of pathological behaviour. By 
contrast particle trajectories in L(3) are well behaved and this approximation gives 
results in excellent agreement with the exact top-hat solution as long as the Void is 
not too underdense. For very underdense Voids, L(3) evacuates the Void much too 
rapidly leading us to conclude that the Zel'dovich approximation L(l), remains the 
best approximation to apply to the late time study of Voids. The behavior of high 
order approximations in spherical voids is typical for asymptotic series and may be 
generic for Lagrangian perturbation theory. 
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1 INTRODUCTION 

It is widely felt that gravitationally bound systems such as 
galaxies and groups and clusters of galaxies arose due to 
the gravitational amplification of small density fluctuations 
already existing at the time of last scattering (and reflected 
in the angular fluctuations in the CMBR measured by the 
COBE satellite). 

The growth of density perturbations in a homogeneous 
and isotropic Universe can be studied in essentially two dis- 
tinct ways. In the Eulerian approach (pioneered by Jeans 
and Lifshitz), solutions to the Euler-Poisson system are 
sought by means of a dimensionless density perturbation 
5(x, t) — [p(x, t) — po]/po, which is expanded in a perturba- 
tion series having the form 



x = q + *(t,q) 



(2) 



£<5 ( ' l) (x,t) 



£ !>?(*) A<»>(x) 



(1) 



in a spatially flat matter dominated Universe (similar ex- 
pansions are made for the velocity field and gravitational 
potential). corresponds to the linear theory limit for 
which D+(t) oc a(t) oc t 2 ^ 3 . Eulerian perturbation theory 
can be successfully applied as long as the series (ji]) con- 
verges i.e. \S\ << 1. 

An alternate viewpoint to Eulerian theory was proposed 
by Zel'dovich who suggested perturbing particle trajectories 
as a means to understand gravitational instability. In a gen- 
eral Lagrangian framework 



where x and q are the Eulerian and Lagrangian coordi- 
nates of a particle/fluid element. The displacement vector 
\E' is related to the initial velocity field and its derivatives, 
in Zeldovich's approximation = = D+(t)tp^ (q), 

ip^ being the initial velocity field. In recent years attempts 
have been made to consider higher order corrections to the 
Zel'dovich approximation by treating "J/ as an expansion 
(Moutarde et al. 1991, Bouchet et al. 1992, Buchert 1993, 
Catelan 1995, Bouchet et al. 1995) 



i = l 



(n) 



Expansions of the form (Q), (^) are frequently referred to 
as n-th order Lagrangian perturbation series L(n), the first 
order series L(l) being the Zel'dovich approximation. One 
might expect (in analogy with Eulerian perturbation the- 
ory - E(n)) L(n) to get more accurate with increasing 'n' as 
higher order terms are included in *P (Coles & Sahni 1995). 
This is indeed the case during the weakly non-linear regime 
as demonstrated by Munshi, Sahni & Starobinsky (1994) 
who showed that results derived from L(n) matched those 
derived from E(n) as long as the system is weakly nonlinear. 
An important issue which has so far remained largely un- 
addressed concerns the domain of convergence of the L(n) 
series (^) & (^) . Although several recent studies have demon- 
strated the increasing accuracy of higher order L(n), most 
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treatments (with the exception of Bouchet et al. 1995) simu- 
lated overdense regions and left unanswered the related (and 
important) issue of the accuracy of L(n) in underdense re- 
gions (Voids). We shall address this issue in this letter by 
studying the behavior of L(n) in spherical underdense re- 
gions complementing and extending the analysis of Munshi 
et al. (1994) for overdense regions. (We leave a more general 
analysis treating generic initial conditions to a later work 
(Buchert, Sahni & Shandarin 1995)). 



2 LAGRANGIAN PERTURBATION THEORY 
IN VOIDS 

(n) 

In a spatially flat, matter dominated Universe factorise 



D+(t) = ^^a(t)^t 2/s 



(4) 



Furthermore higher orders in vj/'™' (n > 1) can be con- 
structed from lower orders by means of an iterative pro- 
cedure (Moutarde et al. 1991, Buchert 1992, Lachieze-Rey 
1993), so that Vf 5 = -<9<Mq)/%, for L (l) (the Zel'dovich 
approximation), for L(2) (equivalently the post-Zel'dovich 
approximation) (Bouchet et al. 1992, Buchert 1992) 

3 



^ = -s((0-«) 



^ (2) = 

For L(3): 



(3) 



(5) 
(6) 

(<) 



Here coma implies partial derivative with respect to q, sum- 
mation over repeated indices is assumed (cf. Juszkiewicz et 
al. 1993, Bernardeau 1993). 

Let us consider the case of a spherical top-hat under- 
density for which 5^(x,t) = — a(t), the initial gravita- 
tional potential can be obtained from the Poisson equation 
and has the time independent form (Munshi et al. 1994) 
<f>o = a 2 r 2 6/9t 2 = -a 3 rg/9t 2 . 

The exact solution for the expansion of a top-hat Void 
can be parametrised as 



R(9) 
a(6) 



(3/lO)(cosh0 - 



(3/5) 



(3/4)(sinh0- 9) 



2/3 



(8) 



Here R = a(i)r/ro is the physical (and r the comoving) 
particle trajectory. Knowing r(a) the equation of mass con- 
servation allows one to determine the density contrast 

5 = (r /rf - 1 (9) 

which for negative density fluctuations has the form 



. ... 9(0-sinh0) 2 
Oth [0) = - 



(10) 



2 (cosh 9 - l) 3 

The related dimensionless peculiar velocity (1/aH) div x u = 
9 (u = ax, x is the comoving coordinate) may be determined 
from the Eulerian mass conservation equation (Munshi et al. 
1994) 



a-^-5 + (1 + 5)9 = 

da 



(11) 



Particle trajectories for spherical top-hat expansion are 
easy to derive in L(n), substituting 0o = -a 3 r 2 /9t 2 in (§) 
& (0) we get 



ri (a) = r (l + -) 



r 2 (a) = r 1 + 



3 21 



, a a 2 23a 3 

r 3 (a) =r H 1 

w 1 3 21 1701 



(12) 



r n is the particle trajectory in L(n). From ([[§) and the con- 
servation condition (|9j) we obtain expressions for the density 
in a top-hat void 



i a a , 
' ' • 3 21, 1 ' 

, a a 2 23a 3 \" 3 , 
53=11+ 3-2T + 170t) - 1 



(13) 



The related expressions for the velocity field are determined 
from dill) 



(a- 



n = a [ 1 + - 



2 2W,, , a a 

r )/(1+ 3-2i } 



7^ + 189 a)/(1+ 3 



21 + 1701 



23 3s 
a ) 



(14) 



(The expressions for a top-hat overdensity can be ob- 
tained from (^) - ( |l4| ) by the transformation a — > —a.) 

Our results for the evolution of overdense and under- 
dense regions are summarised in figures (|l| & ([]) for top-hat 
density and velocity fields respectively. Looking at the right- 
hand side of figure ([!]) (and the left hand side of figure (||)), 
we see that higher order L(n) outperform lower orders in 
matching the results of exact top-hat collapse. The remain- 
ing two panels - showing top-hat Void expansion, however 
give a strikingly different picture. We find that although 
L(2) and L(3) do paint an accurate picture at early times, 
their accuracy declines once the top-hat density falls below 
5th ~ —0.7. In fact, for L(2) particle trajectories begin to 
converge (rather than diverge) at late times, leading to the 
'collapse' of a top-hat Void, a rather pathalogical result ! 
(also see Bouchet et al. 1995). Third order fares consider- 
ably better than second order since particle trajectories al- 
ways diverge leading to Void expansion. However, as we see 
from figures (^) and (|^), the third order Void expands too 
rapidly at late times, causing the asymptotic limit 5 —* — 1 to 
be reached much too early. (Interestingly 52(a) > 5th (a), 
whereas 6i(a), 63(a) < <5th(o).) (In simulations with ran- 
dom initial conditions one might expect shell crossing to 
occur before the density contrast in Voids has dropped to 
5th ~ —0.7, as a result the pathological behaviour of L(2) 
may not be discernable in simulations which evolve parti- 
cles till the epoch of shell crossing but no further (Buchert, 
Melott & Weiss 1994).) 

Although the above discussion has been limited to 
homogeneous Voids, our main conclusions appear to be 
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more general. To demonstrate this we consider the case 
of a spherical density perturbation with an arbitrary ini- 
tial potential ^o(ro). Expressions for L(l), L(2) and L(3) 
turn out to be remarkably simple: tp^ 



'(i) 



-3/7(^/r ) 2 r , $ 



'(3) 



-(<P'a/ r °) ro, 
-23/63(^/r ) 3 f , 



where <fia — dc/>o/dro- As a result 



rl = r o 1 - (D+<j>' /ro) 



r2 = r 



(D +< t>' /r ) - -(D+^/ro) 2 



J"3 = J" 



l-(D + ^/ro)- 7 (D+^/r ) 2 



23 
63 



(15) 



(Related expressions for the density and velocity fields can 
be easily determined from ([]) and (|ll|).) Specialising to 
power law potentials cj>o = —Arg, A — a 3 /9t 2 = constant, 
we get 



T2 = r 



1 + nD+Arr 2 - ^{nD + Arr 2 f 



(16) 



(D+A — a/6). We find that for L(2) the displacement vec- 
tor's i/jl 2 ' and if)!- 1 ' have opposite directions leading to ex- 
pansion at early times being replaced by collapse at later 
epochs. The form of (|E|), ^ suggests that this behavior may 
be generic. Indeed the transformation </>o — > ~<t>o (trans- 
forming a concave potential into a convex one) results in 



t/^ 1 ' — > —ipY', ipY' — > tpY', the fact that the displacement 
field ipl remains inwardly directed regardless of the sign 
of the potential, results in particles flowing inwards at late 
times leading to the eventual 'collapse of Voids' in L(2). 
Thus we find that although L(2) and L(3) start out be- 
ing more accurate than L(l) at the commencement of ex- 
pansion, their accuracy declines with time, and somewhat 
surprisingly, L(l) provides a more accurate picture of Void 
expansion than either L(2) or L(3) at very late times (see 
figures 

It would be interesting to determine whether the para- 
logical behaviour of L(2) afflicts all even orders of the La- 
grange perturbation series or whether it is restricted to sec- 
ond order alone. Although a completely general treatment 
of this kind lies beyond the scope of the present letter, it 
is possible to assess what happens at any arbitrary order 
for the case of a top-hat Void. This case is exactly solvable 
since an expression for the particle trajectory in L(n) can be 
found by expanding the exact top-hat solution (^) in powers 
of a (Munshi et al. 1994), leading to 



= r 



n 



(17) 



(a > 0; for n = 5, ai = 1/3, a 2 = 1/21, a 3 = 23/1701, 
q 4 = 1894/392931, a 5 = 3293/1702701.) The form of 
( fl7j ) indicates that all even order terms will have negative 
signs leading to the eventual contraction of top-hat Voids 
in L(n=even). Odd orders are better than even orders but 
overestimate Void expansion at late times giving rise to very 
empty Voids as demonstrated in figures (jl|) & (|^) for L(l) 
& L(3). 



3 CONCLUSIONS 

We have discussed the evolution of spherical Voids in the 
framework of Lagrangian Perturbation series L(n) of which 
the Zel'dovich approximation is the first order solution. We 
find that L(n) with n-even overestimate the density in Voids, 
whereas L(n) with n-odd underestimate it. On the whole we 
find that L(n=odd) provide better descriptors of Void ex- 
pansion than L(n=even). The Zel'dovich solution L(l) out- 
performs L(3), L(5),... at late times which is typical for 
asymptotic or semiconvergent series. We have not proved 
it but speculate that this may be a generic property of the 
Lagrangian and probably Eulerian perturbative theories of 
gravitational instability. 

It should be mentioned that the relative accuracy of 
L(l) in spherical top-hat Voids has earlier been tested 
against the accuracy of several Eulerian approximations 
(EA) by Sahni & Coles (1995) and Bouchet et al. (1995), 
who demonstrated that L(l) performed significantly better 
than EA at all times. The results of this letter strengthen 
that conclusion and show that the Zel'dovich approxima- 
tion is the best non-linear approximation to apply to the late 
time study of spherical Voids (also see Sahni, Sathyaprakash 
& Shandarin 1994). 
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Figure 1. The density contrast &app m Lagrangian perturbation 
series L(n) is shown plotted against the exact top-hat solution 
&EX for underdense regions (lower left) and overdense regions 
(upper right). We find that whereas the accuracy of L(n) increases 
with 'n' when describing the behaviour of overdense regions, L(n) 
with n > 1 , do not fare as well when applied to underdense regions 
or 'Voids'. For Voids we notice that although L(n) with n = 2, 3 
are initially more accurate than L(l) (Zel'dovich approximation), 
their accuracy worsens with time. We also find that L(2) shows 
pathalogical behaviour at late times when Sex < —0.7. (Note: 
Different scaling has been used to label the negative and positive 
axes.) 



Figure 2. The dimcnsionlcss velocity field 8app in Lagrangian 
perturbation series L(n) is shown plotted against the exact so- 
lution 9ex, for overdense regions (lower left) and underdense 
regions (upper right). We find that L(n), n = 2,3 give better 
results than L(l) for overdense but not for underdense regions. 
(Note: Different scaling has been used to label the negative and 
positive axes.) 
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